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Abstract 

We develop the concept of trajectories in anyon spectra, i.e., the continuous 
dependence of energy levels on the kinetic angular momentum. It provides a more 
economical and unified description, since each trajectory contains an infinite 
number of points corresponding to the same statistics. For a system of non¬ 
interacting anyons in a harmonic potential, each trajectory consists of two infinite 
straight line segments, in general connected by a nonlinear piece. We give the 
systematics of the three-anyon trajectories. The trajectories in general cross each 
other at the bosonic/fermionic points. We use the (semi-empirical) rule that all 
such crossings are true crossings, i.e. the order of the trajectories with respect to 
energy is opposite to the left and to the right of a crossing. 
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1 The concept of trajectories 

Anyons 0, 1 are two-dimensional particles whose wave function ^ obeys the inter¬ 
change conditions 

Pmn^ = exp(ivrz/)^, (1) 

where P mn denotes continuous anticlockwise interchange of particles m and n, such 
that no other particles are encircled, and v is the statistics parameter, which may 
be any real number. For a given Hamiltonian, not explicitly dependent on v and 
with a discrete spectrum for any u, the energy eigenvalues E k (v) and (usually) the 
corresponding eigenfunctions , 0fc( z/ ) will be continuous functions of v. A continuous 
function E k (u), for a fixed index k and v G (— 00 , 00 ), is what we will refer to as an 
(anyon) trajectory. For an earlier discussion of the v dependence of the states, see [3]. 

Since the exponent in (|l|) is the only place where v appears in this formulation, the 
spectrum and the set of energy eigenstates are obviously periodic in v with period 2; 
thus, all the information is contained in the sets {E k (u)} and {ijjkiy)} for v G [0,2). 
However, individual trajectories are not generally periodic in w, one way to see this is to 
observe that when v increases continuously, the kinetic angular momentum L changes 
according to the formula 

m = m + 1 ^(n - i)v, ( 2 ) 

with L(0) an integer]], so that changing v by 2 will not bring us back to the same 
state. It is assumed that the Hamiltonian is rotationally invariant, so that L is a good 
quantum number. Thus, E k (u + 2 £), where i = 0, ±1, ±2,..., will be the energy of 
some state with statistics parameter v, but the index of that state will be different 
from k: 

E k (v + 2£) = E k ,(v)\ (3) 

where k' depends on k and £: k! = k'(k,£). At any given u, the trajectories k and 
k' look different, having different angular momenta L k (v) and L k '(v) = L k (u + 2 E) = 
L k (v) + N(N — 1)£, but when viewed on the interval (— 00 , 00 ), they are seen to be 
copies of one and the same trajectory, shifted by 2 i along the v axis. 

To represent the trajectories in a way which is not redundant, it is convenient to 
take L as a parameter instead of v. L is directly related to v by equation (H), it is 
an observable (gauge invariant) quantity, and the dependence E(L) is the same for all 
the trajectories that are shifted copies of one another. We arrive at the viewpoint that 
the states of an iV-anyon problem may be obtained from the set of trajectories, i.e., 
functions E(L) for L G (— 00 , 00 ). This provides a more economical and unified way 
of decribing the TV-anyon spectra. For each trajectory there is an integer value L(0) 
such that the statistics is bosonic at angular momenta L = L( 0) +N(N — 1)1. We may 

5 This formula can be proved by noting that a 27r rotation of the whole system multiplies the 
wave function by a phase factor, which equals, on the one hand, exp[2i7rL(y)], and on the other hand, 
exp[iV(7V — l)i7w], since each of the N(N — l)/2 pairs of anyons is interchanged twice. 
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always choose —N(N — l)/2 < L( 0) < N(N — l)/2. The value of L(0) groups the set 
of trajectories into N(N — 1) classes. Only trajectories from the same class can cross, 
i.e. have the same E and L at the same statistics. 

A part of a trajectory between any two neighboring bosonic points corresponds to 
what is usually referred to as one anyon state on the interval v G [0,2). It is useful 
to have a geometric picture of this: If one plots a trajectory for — oo < L < oo and 
then wraps the plot around a cylinder of circumference N(N — 1), then points which 
correspond to v differing by an even number (that is, to the same statistics) will fall 
on the same vertical line. On the surface of the cylinder one will see the set of all 
iV-anyon states corresponding to the trajectory. Thus, an infinite number of pieces of 
trajectories (or states in the usual terminology) on the interval v G [0, 2) are shown to 
make up one single trajectory on the interval v G (— 00 , 00 ). In other words, certain 
anyon (and in particular, boson or fermion) states that are usually considered entirely 
different, are in fact parts of one and the same continuous pattern. This is reminiscent 
of the method of organizing the spectra of particle physics into Regge trajectories. In 
fact, by applying the concept of Regge trajectories to a two-dimensional, two-particle 
system, one obtains the same grouping of bosonic/fermionic states. 

It follows that if some trajectory E(L) is analytic, then it is in principle sufficient 
to fold only one of its pieces. The rest is uniquely determined by analytic continuation. 
This means that it is sufficient to fold one anyonic state in some range of v and then an 
infinite number of others are obtained “automatically” by analytic continuation and 
periodicity. However, trajectories are not always analytic. With our assumption that 
the Hamiltonian does not explicitly depend on v (and involves only non-singular inter¬ 
actions) both energies and wave functions will depend analytically on v, except possibly 
at bosonic points. The source of non-analyticity is that the relative angular momentum 
of one pair of particles becomes zero, which can only happen at bosonic points. When 
there is non-analytic behavior at some point with degeneracy, one must determine how 
each trajectory continues through this point by investigating the corresponding wave 
functions. They should change continuously across these points. 

2 Two- and three-anyon trajectories 

We are going now to apply this reasoning to the problems of two and three anyons in 
a rotation symmetric harmonic potential. We scale the variables so that h = 1, the 
mass rri — 1 and the angular frequency uj — 1. Two more preliminary remarks are in 
order. First, the potential is parity invariant; hence, if 0 is an eigenfunction of the 
Hamiltonian and satisfies (]T|), then its complex conjugate 0 is an eigenfunction of the 
Hamiltonian with the same energy and satisfies (|T) with —v instead of v. In particular, 
this implies that all information about the states is in fact contained in the interval 
v G [0,1]. 

Second, there is the tower structure of the spectrum [[|, 0 • It has been observed that, 
for any number of anyons, all the states come in towers, the angular momentum being 
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the same for all members of a tower and the energies being E(v), E(u)+ 2, E{v) +4,.. 
where E{y) is the energy of the lowest, “bottom” state. This sequence of levels is due to 
radial excitations. The radial coordinate r is defined by r 2 = {%i — A") 2 + (r/j — Y') 2 , 

with (A", Y) the center-of-mass coordinates. Thus, it is sufficient to find only the 
bottom states (consequently, bottom trajectories); from now on, we will always mean 
these, unless otherwise specified. 

Now we will demonstrate that in the two-anyon problem there is only one (bottom) 
trajectory, namely the continuation of the ground state. Recall that the one-particle 
spectrum in the harmonic potential consists of states with E = 1,2,3,... and L = 
E — 1, E — 3,..., — (E — 3), — (E — 1) (so the degeneracy equals the energy). Define 
complex coordinates by Zj = {x 3 + iyj)/\/2. Since the center-of-mass motion is trivial, 
we will always concentrate on the relative motion only. For two anyons, the complete 
set of (not normalized) solutions may be written as 

V»in(z, z) = z\ 2l+v \ iFi (-n, |2 1 + u\ + 1; zz) exp (-zz/2) (4) 

with z = z\ — Z 2 , l an integer, n a non-negative integer, and 5 standing for z if 21+ v > 0 
and for z if 21 + v < 0; the energy and the angular momentum are 

Ei n [y) = \21 + v\ + 2 n + 1, Li n (y) = Li(v) = 21 + u, (5) 

respectively. Towers consist of states with the same l and different n\ in accordance 
with the aforesaid, it is sufficient to consider the bottom states only, for which n — 0 
and 

E = \L\ + 1 (6) 

irrespective of/. Indeed, all the states in question belong to one and the same trajectory, 
because the relation (|3|) does hold, in the form E 0 (v + 21) = E t (u) (and the same for 
the wave function). For — 1 < v < 1 this is the ground state, 

V-oOM) = exp (-2^/2), (7) 

Eo(u) = \v\ + 1. (8) 

Thus, all bottom states are continuations of one ground state. Equation (g) exemplify 
our previous remark that a trajectory may become non-analytic when the relative 
angular momentum of one pair of particles (here L itself) becomes zero. By adding 
the pair potential g 2 /\z\ 2 to the Hamiltonian, (|J) is changed to E(L) = sJL 1 + g 2 + 1, 
which demonstrates that singular interactions may lead to different behaviours. 

We go now to the problem of three anyons, which is extremely interesting due to 
its nontriviality, on the one hand, and the possibility of a more or less exact analysis, 
on the other hand. Recall some results available. A state starting with (. E , L) at 
v = n (n integer) may reach one of the following points at v = n + 1: (E + 3, L + 3), 
{E + 1, L + 3), (E — 1, L + 3), (E — 3, L + 3) || [jj. In other words, in any interval 
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from a bosonic to a fermionic point or from a fermionic to a bosonic point, each state 
is characterized by an (average) slope 


A E A E 

s = —— = 3 — = ±1 


Av 


A L 


or 


±3, 


( 9 ) 


where A E is the change in energy, and A L = 3Az/ = 3 is the change in angular 
momentum. For the s = ±3 states the dependence E(L) is linear, so that s = dE/dv 
is the slope at any point in the interval, and their wave functions may be written down 
exactly 110- For the s = ±1 states the dependence is nonlinear, so a slope of ±1 
is indeed only an average slope. Now, there are exact expressions for the multiplicities 
of states with all slopes |]TTJ . If b n (E,L) denotes the number of bottom states in the 


relative motion spectrum that go from (E, L ) to ( E + n, L + 3) as v goes from 2 m 
to 2 m + 1, i.e., from a bosonic to a fermionic point, then there are the asymptotic 
expressions for E, L 1 (terms of order unity being omitted) 


b+\E,L) = (!<£<£), (10) 

b*\E,L) = 2E ~\U~ E \ (_| <£<*;), (11) 

and the exact formulas fe~ 1 (F', L) = b +1 (E — 1 , — L — 5), b~ 3 (E,L) = £> +3 (iA-6, — L — 6). 
In each of the formulas, L has to be within the respective interval specified and the 
equality L = E (mod 2) has to hold, otherwise b n (E,L) = 0. Recall also that for all 
states at any statistics, there is the inequality \L\ < E ||T2|| . 

As one goes along a trajectory, the slope will change at certain values of L, so any 
trajectory possesses a sequence of slopes and points of slope change. We make now the 
following statements. 


(i) For each trajectory there exists an L + such that s = +3 if and only if L > L + . 

(ii) For each trajectory there exists an L_ such that s = —3 if and only if L < L_. 

(iii) The change of slope from/to ±3 can occur at bosonic points only, while the 
change from —1 to +1 or vice versa can occur at both bosonic and fermionic 
points. 


To prove statement (i), note first that the inequality E > L and the fact that dL/dv = 3 
force every trajectory to have s = +3 at least for some large positive values of L. Now, 
the explicit form of the s = +3 states [[h|, ii is such that if one of them exists for 
some v, it does exist with s = +3 for any v > u, that is for any L > L = L(v)\ in other 
words, if s = +3 at some point, then s = +3 everywhere to the right of that point. 
This concludes the proof. A “mirror reversed” reasoning proves statement (ii). 

Statement (iii) is proved upon noticing that the change of slope from/to ±3 always 
implies nonanalyticity in the function E(L), because it is a change between a linear and 
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a nonlinear dependence. Such nonanalyticity can only happen at bosonic points, as 
discussed in Sec. 1. This is also associated with the breakdown of regular perturbation 
theory^ at such points || 

Note that since statements (i) and (ii) mean that each trajectory has s = —3 for 
some values of L and s = +3, for some others, it follows that each trajectory is indeed 
nonanalytic at least at one point. (This is of course equally true for the fV-anyon 
problem, where the extreme slopes, with linear behavior, are ±N(N — l)/2, and each 
trajectory possesses both of them just like in the case at hand; cf. Q.) 

So far we may conclude that the behavior of each trajectory has the following 
features: At large negative L, slope —3, at some point (A_,L_) a change to ±1, for 
L_ < L < L + some sequence of slopes ±1 and finally a change to +3 at (E + , L + ) and 
always +3 at L > L + . Note that there is exactly one trajectory for which = L + , 
that is which has slopes +3 and —3 only: the one which contains the ground state near 
Bose statistics, 

{ (z 12 z 2 3 z 3 i)~ l/ for v < 0, 

( 12 ) 

(z 12 Z23Z 3 i) u for V > 0, 

where Zjk = Zj — Zk, and the overall Gaussian factor is understood. 

We will now prove that 

(iv) Each bottom trajectory may be unambiguously labeled by its (E + , L + ) point. 

Indeed, given an arbitrary bosonic point ( E,L ), it is easy to see that the number of 
bottom trajectories for which (E + ,L + ) = (E,L), equals 

~b +3 {E, L) - b +3 (E - 6, L - 6) < 1, 

where the inequality follows from eq. fllUp. The equality is reached, meaning that a 
trajectory with (E + , L + ) = (. E , L) exists, for each (E, L ) such that ^ < L < E and 
L = E (mod 2). Like the equation itself, this statement is true up to terms of the order 
of unity. Fig. 1 shows the exact picture of the distribution of the (E + ,L + ) points on 
the plane (the bullets and triangles will be explained later). 
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3 The linear parts of the trajectories 


We proceed to show how all the linear parts of the trajectories (and thus all the three- 
anyon linear states) are constructed as excitations of the “ground trajectory”, the one 
which contains the ground state. By virtue of symmetry, it is enough to consider 
s = +3 states only. It is convenient to use the coordinates proposed in [[1(J, which are 


6 The usual way to introduce perturbation theory is by replacing the anyonic boundary conditions 
by a v dependent Aharonov-Bohm-type (perturbative) interaction. 
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the discrete Fourier transform of the complex particle coordinates. The relative motion 
is described by 


u = 


-j= (-1 + + h 2 -s) , v = -j= [z\ + rfz 2 + VZ 3 ) , 


V3 


(13) 


where 77 = exp (2in/3) = (—1 + fv / 3)/2. Introduce the creation and annihilation 
operators 


a u = + flu), 

bu 


4 = 


^/_(w d u ), 

(14) 


(15) 


and, with u —> v, the same relations for a v , aj, b v , b' v . Then the relative Hamiltonian 
and angular momentum operator become 


H = ala u + c4a„ + + b\b v + 2, (16) 

L = a} u a u + a*a„ - b ] u b u - b\b v . (17) 

The commutation relations are [a. u ,a,f] = [b u , b\] = [a v ,al\ = \b v ,b\] = 1, and all other 
commutators vanish. Consequently, if -0 is a common eigenstate of H and L with 
quantum numbers (E, L ), then {a\) 1 {b\ l ) m {b\) n ^ is also a common eigenstate with 

(E + k + l + m+n, L + k + l— m — n). However, to yield true anyonic eigenstates, a com¬ 
bination of the creation operators must be fully symmetric and produce wave functions 
that are not singular at the points where the positions of two particles coincide. 

The pairs (u, v) and (h, u) define two equivalent irreducible representations of the 
permutation group S 3 : for example, the pair (u,v) transforms, under the six possible 
permutations, into (u,v), (■ r]u,rfv ), (rj 2 u,r]v), (v,u), (r]v,r] 2 u), (r} 2 v,r}u), respectively. 
Therefore fully symmetric quantities are those of the form u k u l v m v n + v k v l u m u n with 
k — l — m + n = 0 (mod 3). The pairs (a^, a£), {b\ n 6 ],), and (v,u) all transform 
like (u,v). Taking successive products of the above-mentioned representations of S 3 
and decomposing them into irreducible representations, it is straightforward to prove 
that all symmetric polynomials in , a£, b^, b\, can be expressed as polynomials in the 
following basic symmetric polynomials: 


(k 

l m n) 




(0 

3 0 0 ) 

J 

c 3,-3 

= (bl) 3 + (bl)\ 


(0 

1 0 1 ) 

J 

c 2,-2 

= b{bl 


(0 

2 1 0 ) 

J 

C 3 ,-l 

= (bl) 2 al + (bt) 

V 

(1 

1 0 0 ) 

c 20 

= albi + albl, 


(2 

0 0 1 ) 

J 

C 31 

= ( a u) 2 bi + (4) 

2 6t 

(1 

0 1 0 ) 

J 

c 22 

= 


(3 

0 0 0 ) 

J 

c 33 

= (4) 3 + (<4) 3 - 
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(19) 


The meaning of the subscripts is that 

[H,cl q ]=p c l q) [L,c^ = qc\ q , 

so c' pq changes the energy by p and the angular momentum by q. It remains to check 
for possible singularities at u = v. All the +3 states have the form, up to the Gaussian 
factor exp (—uu — vv), 


3 \v 


P{Zi, Zi) [(zi - z 2 )(z 2 - z 3 )(z 3 - Zi)] v = P(u,u,v,v)(u 3 - V s ) 


with P a polynomial and v > 0. Therefore the only source of singularity can be 
differentiation of the last factor. This immediately implies that a^, which involves du 
only, is always regular, and so is aj; hence, c\ 2 and C 33 are regular. Now, 

bi(u 3 - v 3 r = -j=u(u 3 - v 3 r - - v 3 y-' 

shows that b' u by itself may be singular, due to the last term. Nevertheless, 

C2o(w 3 — v 3 ) u = regular terms — —u 3 (u 3 — v 3 ) v _1 + —v 3 (u 3 — v 3 ) v ~ l 

2 2 

3 v 

= regular terms-— (u 3 — v 3 ) u 


is regular. Further, 

c 3 i (w 3 — v 3 Y = regular terms — ^ u 2 v 2 (u 3 — v 3 y~ l + ^ v 2 u 2 (u 3 — u 3 ) 1 ' -1 , 

but _ 3 , _ 2 and c 3 are singular. A generic (not normalized) linear state of slope 

s = +3 may then be written as 

i>Umn(u,V;v) = (4,)‘(4l )' (cj 2 ) m ( cj 3 )"-00 ( «, V, v) 

= ( 4 o ) t ( 4 i )'( 4 2 )”( 4)”(« 3 - ( 20 ) 


with 


E = 2k + 3/ + 2 m + 3n + 3u + 2, 
L = l + 2m + 3n + 3v. 


( 21 ) 

( 22 ) 


(For the —3 states, of course, Cp- q would take the place of c^ pq .) This is a classification 
of all linear three-anyon states (cf. (SI- 0 )- However, some of these states are tower 
excitations and some are different parts of the same trajectories. First, the operator c^q 
is nothing but the tower raising operator. It never produces singularities when acting 
on any state, linear or nonlinear (T^j, thus providing the tower structure of the whole 
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spectrum. Hence put k — 0. Second, note that V’oo mn(u,v,v) = 2 m+3n / 2 (uv) m (u 3 + 
v 3 ) n (u 3 — V 3 y e ~ uu ~ vv a nd consequently 

^0,0,m,n+2(w, v) - 4^0,0,m+3,n («, «5 = 8[(« 3 + U 3 ) 2 - 4(uu) 3 ]^00mn(u, V ; is) 

= 8(u 3 - V 3 ) 2 y 00 mn(u, V, is) 

= 8y 00mn (u,v, is + 2); 


in other words, the wave functions (|20|), being linearly independent for 0 < z/ < 2, are 
no longer so for 0 < is < oo, therefore some of them do not lead to new trajectories 
and should be excluded from the trajectory counting. As the last formula shows, it 
is sufficient to restrict to n < 2. Thus, a bottom trajectory may be labeled by three 
numbers /, m — 0,1, 2,... and n — 0,1. This is completely equivalent to the labeling 
by (E + ,L + ), because it follows from (|H1 ) - (|22|) that 


E + = 3 1 + 2m + 3n + 2, 
L + = l + 2m + 3n, 


(23) 

(24) 


since the point of slope change to +3 here is is — 0, and it is straightforward to see 
that for any point ( E +) L + ) there is no more than one set {Imn} such that these two 
equations are satisfied. 

The set of bottom trajectories being two-parametric (plus a “double degeneracy” 
due to n) is due to the fact that a bottom state is identified by three quantum numbers 
(six degrees of freedom, minus two for the center of mass and one for the tower exci¬ 
tations). Two quantum numbers (say, E + and L + ) identify a trajectory, and the third 
one (L) chooses a point on that trajectory. To compare, in the two-anyon problem, 
only one quantum number, L, is enough to identify a bottom state, consequently there 
is only one trajectory. In general, for the TV-anyon problem, the family of trajectories 
will be (2N — 4)-parametric. Let us return once more to Fig. 1; there are two copies 
of one and the same two-dimensional pattern, made of bullets (n = 0) and triangles 
(n — 1), and increasing m or l by 1 means moving within the pattern (3 units to the 
right and 3 units up, or 1 unit to the right and 3 units up, respectively). 

Thus the counting of trajectories is complete, but to find the wave functions, certain 


modifications of eq. (20) are still necessary. First, we need only such functions for which 
C 2 o '0 = 0, where C 20 = a u b u + a v b v is the tower lowering operator, but since C 20 does 
not commute with the c^’s, even for k = 0 some of the functions ( f20| ) will contain an 
admixture of non-bottom states with the same E and L. To correct for this, it turns 
out to be sufficient to replace (<4 ,)' in (|20| ) by another operator, writing 


^lmn(u,V]ls) = Ci(cl 2 ) m (cl 3 ) n lp 0 (u,v;is) 

oc C}(uv) m (u 3 + v 3 ) n (u 3 - v 3 y e - Uu 


where 


Cl = «'[(al ) 3 - (at; 


31/ 


(25) 


(26) 
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and 


Q = a u bl ~ a v b\ (27) 

Clearly, Cj is symmetric, and it is straight- 


18 


is Sen’s supersymmetry operator 
forward to check that c 2 o0; mn = 0 always [c 2 o commutes with Q and gives zero when 
acting on f(u,v)e~ uu ~ vv ]. Also, Ci changes E and L by the same amount as (d^)*, so 
that the formulas (|23|) and (|24]) remain valid. 

Second, a wave function that is continuous at the point ( E + ,L + )—and thus cor¬ 
responds to a trajectory—is, in general, not one of the form (|25|) but a linear combi¬ 
nation of several functions of that form. Indeed, for any fixed l, m, n and 0 < v < 2 
we have (d + 1) degenerate s = +3 states {'ipi : m- 3 j,n(u, v; v + 2 j), j = 0,1,..., d}, 
where d = [m/3]. These are linearly independent, but not orthogonal; and which of 
their linear combinations to choose as a basis in the relevant (d + l)-dimensional sub¬ 
space, makes no difference anywhere except at the point v = 0. Here, the states with 
j = 1,..., d continue as linear to v < 0, and the basis in the relevant d-dimensional 
subspace may still be chosen at will, but one state “decouples”, i.e., becomes nonlinear, 
and that is the one orthogonal to this subspace. In general, it will be not 0) 

but a sum J2j=o a j' l /’i,m- 3 j,n(u, v, v + 2j), where the numbers ay/ao, j = 1,..., d, are 
determined by the d orthogonality conditions. Consider the simplest example l = 0, 
m — 3, n — 0, which corresponds to the bosonic point (8,6), see Fig. 2c. The two 
relevant states of the form (|25|) are 0/ = iJj 03 o(u,v] u) and 0 u = il) 000 (u,v] u + 2). At 
v > 0, it makes no difference which two linear combinations of them to choose; but 
at v = 0, "0// continues to be linear, while the state which actually “decouples” and 
continues to v < 0 as nonlinear, is the one orthogonal to 0// at v — 0. An elementary 
calculation shows that it is 0/ + jj0//- Such reasoning has to be repeated at each 
bosonic point with m > 3, and general expressions for the coefficients aj apparently do 
not exist. 


4 The nonlinear parts of the trajectories 


We have shown that each trajectory is uniquely identified by its (E + ,L + ) point; the 
point being given, one should in principle be able to reconstruct the behavior of the 
whole trajectory, and in particular to find all the bosonic and fermionic points it passes 
through. In fact, it would be sufficient to find (E _, L_) and the behavior between that 
point and (E + , L + ) [where E(L) is nonlinear]. To see whether this can be achieved, we 


plot the low-lying trajectories, using the exact multiplicities given in [TIJ, on Fig. 2a-f. 
(There are six different values of T(0) and consequently six plots, one for each class of 
trajectories.) The nonlinear pieces are shown schematically only. The continuity of the 
trajectories is confirmed by the fact that at bosonic and fermionic points, where they 
cross, the total number of trajectories coming from the left and from the right is always 
the same. The main question is how to identify which piece on the left is a continuation 
of which piece on the right. If at a certain crossing point perturbation theory works 
and the first-order corrections to energy are different for all the states involved, then 
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the rule is that the order of trajectories, by increasing energy, on the right is opposite to 
that on the left (because the first-order corrections have opposite signs). This argument 
does not work (a) if the first-order corrections to some of the states are equal and (b) 
when perturbation theory breaks down, at certain bosonic points, as mentioned above. 
A way out then is to find the wave functions numerically, using the method of [|T6] 
and to identify them by comparing their limits as the crossing point is approached. 
As it appears, the above rule—the trajectory which is the n-th from above on the left 
is the n-th from below on the right—holds for all crossings analyzed. In other words, 
all crossings appear to be true crossings. When this rule is used, it becomes possible 
to identify the trajectories completely without actually finding the wave functions. It 
turns out that for all the trajectories there is at most one point of slope change, apart 
from (E + ,L + ) and ( E_,L _), which we denote by (E 0 ,L 0 ), the change being from —1 
to +1 as L increases through L 0 ; clearly, L_ < L 0 < L + . A generic behavior of a 
trajectory is, therefore, 


-3 


( £ _ iL _)4(B 0 ,L 0 )±i(£ + ,L + )±5 


(numbers on top of the arrows meaning slopes), but in fact it may be L 0 = L_ or 
L 0 = L + (for the ground state, L_ = L 0 = L + = 0). 

Table 1 shows the points of slope change for all the trajectories with E + < 12. 


E 

L 

E 0 

L 0 

E + 

L + 

9 

-5 

7 

1 

9 

7 

10 

-8 

8 

-2 

10 

4 

10 

-6 

8 

0 

10 

6 

10 

-4 

8 

2 

10 

8 

11 

-3 

10 

0 

11 

3 

13 

-7 

10 

2 

11 

5 

13 

-11 

9 

1 

11 

7 

7 

-3 

7 

-3 

11 

9 

14 

-12 

10 

0 

12 

6 

14 

-10 

10 

2 

12 

8 

14 

-8 

10 

4 

12 

10 


E 

L 

E 0 

u 

E + 

L + 

2 

0 

2 

0 

2 

0 

6 

-4 

4 

2 

4 

2 

7 

-5 

5 

1 

5 

1 

5 

-3 

4 

0 

5 

3 

4 

-2 

4 

-2 

6 

4 

11 

-9 

7 

3 

7 

3 

5 

-1 

5 

-1 

7 

5 

8 

-4 

7 

-1 

8 

2 

8 

-2 

7 

1 

8 

4 

8 

-6 

6 

0 

8 

6 

9 

-7 

7 

-1 

9 

5 


Table 1. The behavior of the low-lying trajectories. 


The regularity is not obvious here, but it becomes visible when one looks at sufficiently 
high-lying trajectories. Introduce the quantities 


L+ — L 0 L 0 — L_ 


that is, the “numbers of revolutions”, or the numbers of bosonic points that a trajectory 
passes through having slopes +1 and —1, respectively; these may be integer or half¬ 
integer, as L 0 may correspond to either a bosonic or a fermionic point. The following 
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obvious relations hold 


Eq = E + — 2n + , Lq = L + — 6 n_, 

E_ = E + — 2 (n + — n_), L_ = L + — 6 (n + + n_). 

The results, deduced from the numerical analysis, are the following: 


(29) 

(30) 


E+ T ^ E+ 3 L + — E + E + — 2L + . . , A 

— < L + < — : ra+ = ---, n_ = ---, L 0 is fermionic, (31) 


^ < L + < E + : 
2 + + 


E 4 


E 4 


n, = 


n_ = 


Lei is bosonic. 


(32) 


The formulas, as well as the inequality signs, are valid asymptotically for E 1, in the 
same sense as (|10D—(pTl]). One can verify the compatibility of ( pTT|) (|T2|) with (|T(3|)-([IT|); 
also, a simple consistency check comes from the symmetry requirement—a trajectory 
(. E_,L _) —> ( E 0 ,L 0 ) —> (E + ,L + ) must have its partner (E + , — L + ) —> ( E 0 ,—L 0 ) —> 
(E_,—L_): this is satished as well. Note that n + depends continuously on L + , but 
there is a discontinuity in n_ at the point L + = E + /2. Also, the exact results exhibit 
certain periodic structure with period 8 by E + and 4 by L + . Explaining these features 
remains an open item. 

A semiclassical interpretation—rather rough, although—of some features of the 
trajectories may be adduced. When two particles are close together and the third one 
is far from them, that is, say, |pi 2 | < \ps\ [pn = Zi — z 2 , pz = {z\ + z 2 )/2 - ^ 3 ], 
one gets two independent anyonic oscillators, one (pi 2 ) with the statistics parameter 
z/, the other (p 3 ) with the parameter 2v |5]. ^]. The energy in this approximation is 
|2Z i2 + u\ + 2n 12 + 1 + 12Z 3 + u\ + 2n 3 + 1, where in general l 3 / 12 ; as v increases 
through —2li2, the slope changes from —3 to —1 (if l 3 < 0) or from +1 to +3 (if l 3 > 0); 
the other two points of change, however, can not be described correctly, nor can the 
formulas (0M0 be explained. Perhaps a more accurate semiclassical approximation 
would be able to explain them. 


5 Conclusion 

We have developed the concept of trajectories for anyons, noting that many different 
many-anyon states are in fact continuations of each other. We discuss in particular 
the harmonic oscillator external potential. Apart from the center-of-mass and tower 
excitations, in the two-anyon problem there is only one trajectory. For the three-anyon 
problem, we have worked out the classification and the main features of the trajectories: 
The slope of a trajectory always changes as —3 —> — 1 —> +1 —> +3 (where the +1 
and/or —1 pieces may be missing), and the point (E + ,L + ) of the last change may be 
used to label the trajectory. The wave function corresponding to the linear part of 
the trajectory can be written down exactly by applying excitation operators to the 
ground trajectory. Concerning the nonlinear part, we conjecture the formulas, based 
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on a numerical analysis, expressing the lengths of intervals of slope —1 and +1 in terms 
of (£+,£+)• 

Qualitatively, the picture will be the same for the TV-anyon problem: Again, each 
trajectory will have the points (. E + ,L + ) and such that its E(L) dependence 

is linear with slopes ±N(N — l)/2 for L > L + (L < L_) and nonlinear in between. 
Here the family of trajectories will be ( 2N — 4)-parametric. It is quite plausible that 
in this case as well, the sequence of slopes will be regular, from —N(N — l)/2 up to 
+N(N — l)/2 in steps of two. However, to gain a more precise understanding of the 
behavior of the trajectories remains an open problem. 

S. M. would like to express his sincere thanks to the Institutt for fysikk at the NTH, 
Trondheim, and the theory division of the Institut de Physique Nucleaire, Orsay, where 
parts of this work were done, for the kind hospitality and support. 
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Figure captions 

Fig. 1. The distribution in the ( E,L) plane of the points (E + ,L + ) at which the 
trajectories turn linear with slope s = +3. The bullets define the set E + = 3^ + 2771 + 2, 
L + = £ + 2m, and the triangles define the set E + = 3£ + 2m + 5 , L + = £ + 2m + 3, 
with l — 0, 1 ,... and m = 0, 1 ,... in both cases. 

Fig. 2, a-f. The six classes of trajectories in the (E, L) plane, with L(0) = 
—2, —1, 0,1, 2, 3; solid and dashed vertical lines mean Bose and Fermi statistics, re¬ 
spectively, dashes show the multiplicities of the linear states. 
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Fig. 2a 







Fig. 2b 







Fig. 2c 







Fig. 2d 
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